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Ion irradiation during film growth has a strong impact on structural properties. Linear stabil- 
ity analysis is employed to study surface instabilities during ion-assisted growth of binary alloys. 
An interplay between curvature-dependent ion-driven and deposition-driven instabilities is inves- 
tigated. We demonstrate that ion irradiation of growing binary alloys leads to the formation of 
composition-modulated surface patterns. It is shown that the ion-to-atom arrival ratio R is the 
pattern control parameter. Close to the instability threshold we identify different regimes of insta- 
bilities driven by ion- or deposition-induced surface roughness processes, or roughness-composition 
feedback interactions. In particular, the synergistic effects of the curvature-dependent displacement 
and deposition coupling to the preferential sputtering or to the preferential diffusivity are found 
to induce instabilities and pattern formation. Depending on the film growth and ion-irradiation 
conditions, the instabilities show stationary or oscillating behavior. The latter one is exclusively 
connected with ion irradiation. The corresponding phase diagrams are presented in terms of exper- 
imentally accessible parameters. This shows an alternative way to control surface patterning and 
to grow three-dimensional laterally or vertically ordered nanostructures. 

PACS numbers: 81.16.Rf, 79.20.Rf, 81.15.Jj, 81.15.Aa, 68.35.Dv 



I. INTRODUCTION 

Phase separation occurring during thin film growth re- 
sults in a remarkably large variety of morphologies 
It takes place under 'frozen' bulk conditions^ where the 
bulk diffusion is restricted due to low growth tempera- 
tures. Such kinetic constraints for phase separation pro- 
cesses provide unique opportunities to influence the re- 
sulting structure down to the nanoscale. Interface and 
size effects of nanostructured materials synergistically act 
to influence the properties on the macro-scale. For phase 
separated, or nanocomposite films this leads to new prop- 
erties and advanced (multi)functionality which cannot be 
predicted from the film constituents alone^^— There- 
fore, the control over the film nanostructure is of utmost 
importance. 

Thin film growth includes many interacting kinetic 
processes (surface and bulk diffusion, repeated nu- 
cleation, shadowing, surface reaction, growth rate, 
etc.) which are influenced by external parameters 
such as temperature, growth rate, presence of ener- 
getic species, etc£*^~— Structure zone models have been 
constructed to take into account the effects of temper- 
ature and thin film deposition assisting energetic ion 
bombardment These models predict what char- 

acteristic structures will grow such as columnar grains, 
V-shaped columnar grains or equi-axed grains, etc, 
in certain experimental parameter ranges of thin film 
growth. Of particular importance are the surface insta- 
bilities because they can lead to composition nanopat- 
tern formation. 25 During the growth, such nanopattern is 



constantly buried by randomly depositing species. These 
species are again restructured due to the instability. In 
such a way the 2-dimensional (2D) structure appearing 
on the surface is transferred into the bulk, resulting in 
the formation of a 3-dimensional (3D) ordered hetero- 
geneous structure. Considering the above, such a 3D 
nanostructure offers opportunities to control macro-scale 
properties. Therefore, the possibility to tune such surface 
instabilities is one of the paths to control the structure 
of the resulting films. The film growth might give rise to 
one of the following types of instabilities^ 5 -: 

• dynamic instabilities due to interaction of imping- 
ing species with the film surface) 26 ' 27 

• kinetic instabilities if the growth of the surface is 
faster than its equilibration; 25 

• thermodynamic instabilities where the desired 
structure is thermodynamically unstable^ 5 - 

• geometric instabilities where the surface roughness 
shadows the growth of the material^ 5 . 

The analysis of the decomposition kinetics of phase- 
separating binary systems during the film growth has 
demonstrated that the resulting characteristic periodic- 
ity depends on the interplay of surface diffusivity and 
film growth rate^ i 10 ' 28 " — If the growth rate surpasses a 
critical value any phase separation is suppressed ! 28 ' 29 By 
performing linear and non-linear stability analysis, it has 
been shown that in the case of spontaneous phase separa- 
tion surface roughness can couple to composition i 29 ' 29 i 31 
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This phase separation can generate an elastic field which 
acts back on the film morphology^! On the other hand, 
substrate mismatch and compositionally generated stress 
can also induce decomposition in otherwise thermody- 
namically stable alloys J£r— Kinetic effects such as dif- 
ferences in surface atomic mobilities can have a similar 
result. 38 

In such a context, the possibility to use some external 
easily controllable factors whose influence and strength 
on the surface atomistic processes occurring during film 
growth would be comparable or even greater than the 
ones of 'intrinsic' processes is of particular importance. 
The energies delivered by ions to the (sub)surface atoms 
via atomic collisions far exceed the thermal budget of 
the atoms within the irradiated material. This creates 
unique conditions, which cannot be achieved by any other 
means. Ions provide an energy per material atom in 
the range of ~ 1 — 10 eV/atom or more. That is much 
larger than the energy provided by any other processes 
such as solidification, recrystallization, interface energy 
minimization, thermal activation, phase transformation, 
etc. which are around or below ~ 1 eV/atom^ Also 
the directionality is provided which can be externally 
controlled. 39 Assisting ion irradiation is used to influence 
the structure of thin films.— "1 It is well established that 
ions have a huge impact on structural properties like den- 
sity, grain size or texturej 40 ' 41 The control parameter is 
the ion-to-atom arrival ratio R.— 

It is well established that ion irradiation of surfaces in- 
duces pattern formation. 42 For multicomponent materi- 
als such roughness pattern can couple to the composition 
resulting in compositionally modulated surface ripples^ 3 - 
or nanodots.— During ion erosion, the pattern formation 
depends on the compositio n 43 ' 44 which is determined by 
the initial material composition or co-deposition rate of 
impurities.^ 5 - Therefore, one can expect that such an ion 
irradiation also will affect the surface composition of the 
growing multicomponent thin film and will induce pat- 
terning of the growing surface. To our knowledge there 
are no theoretical studies of the ion effects on the sur- 
face compositional and spatial distribution during film 
growth. Such an approach would present alternative 
ways to grow 

• nanostructured surfaces with tunable surface 
roughness/composition patterns composed from a 
material different than that of the substrate; 

• 3D compositional nanopatterns or nanocompos- 
ites with tunable structural properties (periodicity, 
composition, tilt...). 

In this paper, we employ the linear stability anal- 
ysis on the surface roughness and composition during 
ion-assisted bi-component film growth. From the ion- 
induced effects we consider only those which are related 
to sputtering 4 ^ and ballistic-induced surface redistribu- 
tion fluxes42r— Similarly we consider dynamic effects 
from the deposition j 26 ' 27 In order to highlight the effects 



of ion irradiation we consider only the case of an alloy de- 
position, i.e. there are no thermodynamic factors which 
could result in phase separation . 25 ' 37 

The instability analysis is carried out in 'frozen bulk' 
approximation where majority of changes occur at the 
advancing film surface resulting in the 'frozen' structure 
in the bulk.— In general, surface processes dominate in 
the temperature range T <~ 0.5T m ( T - substrate tem- 
perature, T m - melting point) Thus we focus on the 
growth situations where bulk atomic mobility is negli- 
gible in comparison to that of the surface. This is the 
main condition to transfer 2D patterning of the surface 
into the bulk. Otherwise, if the bulk evolves on itself, 
this patterning is lost. In this way the changes in the mi- 
crostructure across the film cross-section represents the 
historical evolution of the microstructure of the film ji2 

From the thermodynamic forces we will consider only 
those driven by surface concentration and curvature 
gradients, since these two are usually used in the lit- 
erature concerning the pattern formation during ion 
irradiation j 42 ' 43 ' 45 The first counteracts the formation of 
composition differences while the second relaxes the sur- 
face roughness modulations. 

We will not consider ad-atoms as separate surface 
species and no separate rate constant is introduced to 
describe the process of surface ad-atoms losing their mo- 
bility on the surface and their incorporation into the fi- 
nal location^ We will assume that surface atoms are 
available for diffusion as long as they are not covered 
by depositing fluxes. This is particularly valid for ion- 
assisted deposition as any atom can be knocked out of 
stable positions to become a mobile species on the sur- 
face. Therefore the characteristic time for diffusion (and 
any other induced process on the surface) is the interplay 
between the growth rate and surface diffusivity* 1 ^ 

Throughout this paper we intend to explore the insta- 
bilities due to ion irradiation. The control parameters 
which are easily accessible from the experiments are the 
ion incidence angle 9 and ion-to-atom arrival ratio R. 
It will be shown that these parameters drive the sys- 
tem from dynamic equilibrium to instabilities during the 
alloy growth. We will omit geometric effects due to de- 
position and consider the perpendicular depositing atom 
incidence which simplifies the analysis. On one hand the 
angular dependency of deposition is described well in the 
literature (see for example Refs. I25ll53h . On the other 
hand we intend to separate these effects from the ion- 
beam induced ones. We will demonstrate that ions in- 
duce composition-modulated surface-roughness patterns 
which show explicit dependencies on R. This originates 
directly from the form of continuum equations describing 
the ion-assisted alloy growth. The instabilities can be of 
stationary or oscillating type. 

In the next section we develop the model of ion-assisted 
film growth which couples height and surface concentra- 
tion of the binary alloy. In section III, the equations are 
linearized and the solutions are discussed with respect 
to time-growing modes. In section IV, different types of 
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instabilities are identified. Their conditions are linked to 
certain combinations of experimental parameters which 
reduces appreciably the parameter number. Section V 
summarizes the different cases in a mathematical form 
and develops the corresponding phase diagrams. The ex- 
perimental realizations of the pattern formation are then 
discussed. Finally, section VI summarizes the findings of 
the paper and presents the conclusions. 

II. MODEL DEVELOPMENT 
A. Evolution equations 

Let us consider a surface in the (x, y) plane with the lo- 
cal height h(x, y). It exhibits the atomic density ps(x, y) 
and consists of the species A and B with the local surface 
atomic fraction eg y), where i — A,B. The deposit- 
ing flux F = Fa + F B of two species A and B (see Fig. [T]) 
is assumed perpendicular to the film plane. The atomic 
surface ratio of species i is cqj — Fi/F and the arriving 
atoms of species i stick to the surface with the probability 
Si. Further we summarize the atomic flux 

Jat = (FaSa + FbSb) /ps ~ (co,a + co,b) FS/ps 
= FS/ps (1) 

where co,s = 1 — Co,a arL d we assume that the mean 
sticking coefficient is independent of species S ~ Sa ~ 
Sb and of local surface concentration. From now on we 
denote the fraction of the species A in the incoming flux 
as cq. 
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FIG. 1: The schematic picture of the experimental setup of 
depositing atomic fluxes Fa and Fb for two species under a 
bombarding ion flux / creating a sputtering yield Ya and Yb 
for the two species. 

The growing film surface is bombarded with ions. Ion 
beam irradiates the surface with the ion flux / at the 
angle 9 to the surface normal. The projection Icos(6) 
of the ion flux density vector onto the film plane is 
along the x axis. Similarly to the deposition flux ( 
Eq. (fT])) we normalize the ion flux per surface atom, 
be. jion = Icos(O)/ ps- Irradiating ions induce sputtering 
of the species i characterized by a sputtering coefficients 
Yi, i.e. number of sputtered atoms of type i per ion. The 
actual sputtering yield of the species i with the surface 



fraction cs,i is YiCs,i- The total sputtering yield Y is the 
sum of the partial sputtering rates Y — Yacs,a + Y b cs,b- 

The evolution of the film surface height h(x, y) in time 
t follows from the conservation of matter at the solid 
surfac o 37 ' 45 ! 54 . Therefore let us consider an infinitely 
small surface area E with the height h(x,y). In the fol- 
lowing we consider only linear analysis. Therefore, we 
neglect the difference between the local height normal to 
the surface and the projection on the z axis, since this 
difference is quadratic in Wh(x, y) i 26 i 48 The time change 
of this material volume h{x, y)E is considered here due to 
the adsorption of atoms by atomic flux and the sputter- 
ing by ions as described above and the redistribution due 
to surface fluxes. Note that in contrast to the derivation 
presented in Ref. H3 we do not assume a speed of the 
moving surface but obtain the speed from the competi- 
tion of deposition, sputtering and relocation fluxes. 

In the 'frozen' bulk approximation, all the processes 
take place within the 'active' volume — AE given by 
the monolayer thickness or atomic diameter A beneath 
the area E. Therefore the material volume is increasing 
with the atomic adsorption by the number of atoms per 
second ~ fijat with Eq. (JTJ) and decreasing with the ion 
flux and sputtering yield ~ — Qji on Y. The rearrangement 
due to surface currents js is given by the change of the 
number of atoms N in the active volume ~ CIN. The 
local fluctuating atomic density n = N/D, changes as 
h = — V • js/E where the current conservation was used. 
Similarly to the case of ion erosion with co-deposition of 
impurities 4 ^, collecting these contributions and dividing 
by the surface area E with the help of Eq. ([T]) one gets 

— = j at A-j ion AV- A 2 V-j s 
at 

= A FS _ A (l^ly + AV . is ). (2) 

Ps \ Ps J 

The current js = (ja + Jb) describes the redistribu- 
tion of the atoms on the surface due to various atomistic 
processes. These can be fluxes induced by deposition, ion 
irradiation or surface diffusion and are described in more 
detail below. Such fluxes change the local surface atomic 
fraction cs.a ( c S,b) of the element A (B). Multiplying 
the mass conservation for the species Afi^ = —V • Ja/E 
with the active volume AE one obtains the contribution 
cs.a ~ — AV ■ Ja- Since this process originates from the 
active volume we have to make sure that the space was 
available for this process which means we have to mul- 
tiply with the concentration of empty states (1 — c^a) 
where c^a — cs,A{x,y,h{x,y) — A) is the atomic frac- 
tion of the element A at the underlying layer. In the 
same way the effect of sputtering is considered as a ki- 
netic process which decreases the concentration cs,a ~ 
—(l — Cb t A)cs,AjionYA and increases it by the complimen- 
tary process of species B, i.e. cs,a ~ 0,^(1 — cs^iion Yb- 

Similarly to Refs. andl45l. by collecting these terms 
together the mass conservation law for the surface atomic 
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fraction cs,a{x, y) of type A atoms results into 

dcs,A 



at 



ja.t(co—C8,A)-(l — Cb,A) (cs,AjionYA + AV • j A ) 



C h , A [(1 - C S .A)jionY B + AV • j B ] 



(3) 



Please, note that the first term represents the source due 
to atomic deposition?^ The deposition of the atoms of 
the type B (A) on the atoms of the type A (B) reduces 
the local surface atomic fraction cs.a { c S,b) of the ele- 
ment A (B), as the surface atoms of the type A (B) be- 
come the part of the bulk. The surface composition does 
not change if the same type of species are covered by the 
atoms from the vapor phase. On the other hand, sputter- 
ing or surface atomic redistribution removes the atoms of 
the type B (A) at the surface location (x,y). The atoms 
from the underlying layer h(x, y) — A become atoms of 
the surface. The local surface composition changes only if 
the uncovered atoms from the bulk are of the type A (B). 
More detailed derivation of this balance can be found in 
the literature! 43 ' 45 ' 54 Note, that in the presence of sput- 
tering (with or without co-deposition) an altered layer 
forms close to the surface which differs in composition 
from that of the bulk i 45 ' 54 Therefore, the consideration 
of the altered near-surface composition depth profile in 
relation to the bulk concentration is necessary in deriv- 
ing the correct concentration balance! 45 ' 54 Only in the 
particular case when the deposition rate is larger than 
the re-sputtering rate, the concentration distribution as 
a function of depth becomes homogeneous. 54 This allows 
a considerable simplification of the equations and will be 
considered below. 

The first term on the right of Eq. © describes the 
local composition changes due to deposition from the va- 
por phase. The second term describe the decrease rate of 
the atoms A due to sputtering and surface redistribution 
as well as deposition fluxes of the atoms A, respectively. 
The third term describe the increase rate of the atoms A 
due to uncovering from the underlying layer at the depth 
h(x, y) — A resulting from the sputtering and surface re- 
distribution as well as deposition of surface atoms of the 
type B, respectively. 



B. Planar film growth 

For the planar homogeneous film growth and absent 
surface redistribution fluxes js, the surface height h and 
atomic fraction ca are independent of (x, y). The balance 
equations Eq. ([2} and Eq. ^ result in the following 
system of equations 



dh? h 



at 

9c s ,a 

at 



- = j at A - j ion A [c s ,aY a + (1 - c s ,a)Y b ] = V 

= jat(c° - C S ,a) 

- jion [CS,^(1 - Cb,A)YA - (1 - CS,A)Cb,AYB\ (4) 



where the coefficients denote the unperturbed, i.e. flat, 
film growth mode and V is the film growth rate. Remem- 
ber that we have abbreviated for the atomic flux j at = 
(Fa + Fb)S/ps, the projected ion flux j ion = Icos(0)/ps, 
£b,A — cs.A{h — A) is the concentration in the underlying 
layer or 'bulk'. 

For steady state conditions there is a growth of a ho- 
mogeneous film, i.e. h = ho + Vt. The element dis- 
tribution becomes homogeneous over the depth?^ 4 - i.e. 
c A (h — A) = Cg A . Then the second equation of Eq. flU 
yields 



f 



-S,A 



2RY 



1+RY- v 7 ! + RY 2 - 4c RY (5) 



where Y = Ya — Yg > describes the preferential 
sputtering of one specie compared to the other one, 
where without loss of generality we assume Ya > Yb- 
The ion-to-atom arrival ratio R = ji on /jat = (FaSa + 
FbS b )/Icos(9) is already visible as the crucial order- 
ing parameter. This stationary solution is valid only for 
V > 0, i.e. when the growth rate is larger than the 
re-sputtering rate. From now on, we will denote that 
the surface concentration refers to the element A, i.e. 



„o _ „o 

C S — C S,A- 



C. Deposition-induced surface redistribution fluxes 

The surface height h(x,y) evolution due to deposition 
in Eq. ([2]) can be described by currents j dep represent- 
ing the dynamic interactions of incoming film forming 
species with the film surface and non-thermal particle 
rearrangements 5 - 5 -^ after the impact £tL Therefore the de- 
position fluxes Fi, the sticking coefficient S, the surface 
density ps and the current terms on the right hand side 
of Eq. © can depend on the surface slope and the 
curvature i 26 ' 27 Let us develop their contributions up to 
second order derivatives of the surface height. 

Particles in the beam moving perpendicular to the sub- 
strate can be deflected when they come very close to 
the surface. This is because the atom-substrate inter- 
atomic forces act perpendicular to the surface, not the 
substrate*^ More atoms arrive at the surface places with 
V 2 /i < that at V 2 /i > 0. This results into the so-called 
'surface steering' effect^ 7 - meaning that the deposition 
flux is not homogeneously distributed over the surface? 5 - 7 - 
The effect is more pronounced for lower atom kinetic en- 
ergies and grazing incidences, but can also be active even 
for normal incidences? 5 ^ Leaving only linear lowest-order 
terms, the deposition flux Fi for the species i can be 
expanded 2 ^ 



^ =F?(l + asT,iV 2 h(x,y) + ...) 



(6) 



where the constant asr.i < describes the dependence 
of the deposition flux Fi on the surface curvature. Such 
a dependence results in a larger depositing atom flux at 
the crests than at the depressions. This induces a topo- 
graphic instability. 
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The sticking coefficient S can also depend on the lo- 
cal curvature. For example, for fcc(lOO) or bcc(lOO) ge- 
ometries with four-fold hollow (4fh) adsorption sites, the 
density of 4fh2£ is expected to be higher at local minima 
V 2 h > than at local maxima V 2 h < 0££. One can argue 
that this effect is not significant for amorphized surfaces, 
as they can easily relax to a random structure indepen- 
dent of the slope and curvature. Considering only linear 
lowest-order terms, this yields^ 7 - 



S = S° (l + a s V 2 h(x,y) 



(7) 



where ag > describes the dependence of the sticking 
coefficient S on the surface curvature. Such a dependence 
of the sticking coefficient on the surface curvature has a 
stabilizing effect on the surface topography. 

If the adsorbing atoms do not hit directly an adsorp- 
tion site, they still posses an ad-transient mobility with 
the significant component in the direction of the initial 
impingement and are 'funneled' downwards until they 
reach such an adsorption site^ Downward funneling is a 
non-thermal process. It induces a downward flux also de- 
scribed by the second term on the right hand side of Eq. 
@. Leaving the linear lowest-order terms one obtains 5 - 6 - 



'n a DFi Vh(x,y) 
Ps 



(8) 



where aoF,i > describes the mobility of funneled atoms 
and pg describes the atomic density of the defect free 
surface. This effect has again a stabilizing behavior on 
the surface topography. 

Deposited atoms can be trapped on the sides of the 
slopes instead of funneling down to adsorption sites. 60 
Such 'restricted downward funneling' can be represented 
by expanding ps^ 



1 

Ps 



P°s 



(l + a SD V 2 h{x,y) 



(9) 



The term clsd < describes the defects which trap the 
funneling atoms and which are dependent on the curva- 
ture. It has a destabilizing effect on the surface topogra- 
phy 

Combining Eqs. ©, ©, ©, © into Eq. © and 
leaving only the linear order terms, one gets 

dh/dt » j at A - j ion A [c s ,aY a + (1 - c s ,a)Y b ] 

-A 2 (V.jl cp + V.j£ p ) (10) 



with 



• dep 



F° S° 

— 1 {asT,i + a s + cidf^ + ci S d) Vft. 

Ps.o 

-CQ.ijutSf^Vh 



(11) 



where Sf op combines the dynamic surface slope and cur- 
vature effects of species i = A,B. Its sign depends on 
the relative magnitude of each of the processes and can 



have either positive (stabilizing) or negative (destabiliz- 
ing) effect. For a given material system, there are only 
few opportunities to control these dynamic effects includ- 
ing the incidence angle. Since we assume the depositing 
flux incidence angle fixed, S dcp will be considered as a 
constant further on. In the following, we will call the 
parameter S dep the curvature-dependent deposition co- 
efficient. 



D. Ion irradiation - induced fluxes 

When penetrating into the solid, energetic ions create 
recoils which retain to the large extent the initial mo- 
mentum of the incoming ions before being stoppedj 47 ' 48 
If the ion beam is oblique to the substrate, then there is a 
component along the sample surface. Due to the incom- 
pressibility of the solid, the excess of the material density 
relaxes into the surfacej 47 ' 48 For the ion flux / with the 
incidence angle 6 from the surface normal, this generates 
ion-induced ballistic relocation currents proportional to 
/cos(6')sin(6'), where the product with cos(#) reflects the 
projection of the ion beam on the surface while the prod- 
uct with sin(6*) reflects the projection of the recoil atom 
flux along the surface. Expansion to linear order yields 
the following expression for the ion-induced relocation 
currents 4 ^— 



;rel 



dh 
dy 



(12) 



where = ftdi is the relocation yield of the species 
i, fi is the number of relocated atoms and di is the av- 
erage relocation distance. The dependence on the ion- 
incidence angle is explicitly included in 5^(0). It can 
change the sign and thus can act to stabilize or destabi- 
lize the surface roughness Assuming a simple angu- 
lar dependence ~ cos(2#) from the model of Carter and 
Vyshniakov, 47 one can see that S™*(0) becomes negative 
for 6 > 45°. This means that the areas with positive 
curvature (d 2 h/dx,y 2 > 0) or depressions start to loose 
the material. The latter starts to accumulate at the ar- 
eas with negative curvature d 2 h/dx,y 2 < 0) or crests, 
i.e. the instability grows. For silicon, it has been demon- 
strated that such ion-induced relocation currents have 
the dominating effect for the ion-induced surface topo- 
graphic instabilities! 49 ' 50 

Part of the ion-displaced atoms gain enough energy 
and leave the solid, i.e. are sputtered. The sputtering 
can depend on the local surface curvature^ To first- 
order this yield s 43 ' 46 



Yi = Y? 



dh . d 2 h 



d 2 h 



(13) 



If the sign of either S^P y (8) becomes negative for a certain 
angular range of 9, it results in larger sputtering yields 
at the depressions than on the crests, i.e. the surface 
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roughness becomes unstable^ While for Si this effect 
has been shown to be weak ) 49 ' 50 the relative strength of 
the curvature-dependent sputtering in relation to the ion- 
induced relocation remains to be determined for other 
material systems. Therefore, for the sake of completeness 
it will be considered in the further analysis. 



E. Diffusive currents 

Surface diffusion fluxes are proportional to the gradient 
of the surface chemical potential fis an d the surface atom 
diffusivity 



diff _ ApsDj 

KrT 



(14) 



Considering only the effects of concentration gradients 
and surface 'stiffness' or capillarity, this leads to 37 ' 43 

ji = -ADipsVcs.i + c Si — — — — V(V) h (15) 

k B T 

where Q, 7, T, and k B stand for the atomic volume, 
surface energy, temperature and Boltzmann's constant, 
respectively. 



F. Final set of equations 

Taking into account Eqs. pU)) - ([T5|) . mass conservation 
(Eqs. © and ([3])) leads to the following set 

= Jat A - j ion A [c s , A Y A + (1 - c StA )Y B ] 
- A 2 V-{j A +j B } 

' '' S ' A = m(c°-Cs,a)-(1-C° s ) \jionCS,AY A + AV • j a] 



Of 



+ c° s [j ion (l - c s , A )Y B + AV • j B ] 



(16) 



where ji = ji dlff + ji rcl + ji Qcp is the sum of the diffusive 
flux, Eq. (|T5|) . relocation, Eq. (TT!?j) . and deposition Eq. 
(jTTJ) . In deriving Eq. (TTBl we used again the fact that for 
film-growing conditions the bulk concentration is equal 
to the steady-state surface value c^a = c s~ 



rcl 1 • dcp 



III. LINEAR ANALYSIS 

A. Linearized equation and parameters 

Now we are going to study the stability of this solution 
to perturbations u(x, y) and <f>(x, y) in height and surface 
concentration, respectively, 

h = h + Vt + u{x, y), cs.a = 4 + <t>& V)- (17) 

Introducing Eq. (jTTJ) into Eq. (fll))) . changing the vari- 
ables r = ja,tt and expressing the height and the spatial 



coordinates in units of the monolayer thickness (U = 
it/A, x — > Ax) we get the following system of equations 



dU 
dr 



8t 



dU 
-Rv— 
ox 

-RY(/) + 



£>V 2 < 



d 2 U 



d 2 U 



dx 2 ' " v dy 2 
-BV 4 U 



RY* 



Rv* 



dx 



R 



. d 2 u 

' dx 2 



s 



-D*V 2 d>- B*V 4 U 



, d 2 u 

V Q y 2 

(18) 



with R = jion/iat and where the coefficients are sum- 
marized in Table HI The seemingly cumbersome combi- 
nation of parameters bear a clear physical meaning. The 
* parameters denote the coefficients in the concentration 
equation of Eq. (ITS)) while all the others are related to 
the equation describing the surface roughness evolution. 

The terms vdU/dx and v*dU/dx in Eq. (fl8]l result 
in a surface drift. As the surface is constantly buried 
by incoming species, these drift terms determine the 3D- 
nanopattern tilt in relation to the sample surface. One 
can absorb the surface drift term ~ dU/dx by a running 
frame transformation x — > x — tRv in the surface rough- 
ness equation of Eq. (TT5|) . However, the drift term in the 
concentration equation will induce nontrivial patterns. 
In the following we neglect the drift terms restricting to 
small drift velocities and devote this analysis to a later 
investigation. 

To illustrate the usefulness of these parameters con- 
sider the diffusion terms D, D* and the capillarity B, B* 
terms in Eq. (TT5|) only. While the D,B* terms can 
change sign due to the difference between species pa- 
rameters and can therefore act to destabilize the system, 
the D* , B parameters are always positive due to sum of 
parameters, see Table fl] In the absence of the deposition 
flux (j at = 0) and ion-irradiation (R = 0), one has always 
D* > D and B > B* and any infinitesimal composition 
or surface roughness fluctuation is suppressed. 

In the following, the term 'displacement' will comprise 
both sputtering and relocations*^ The driving terms 
for instabilities are S x . y and S* y for surface roughness 
and composition, respectively. These terms describe 
the curvature-dependent total and preferential displace- 
ment/deposition rates, respectively. These terms are 
multiplied by the surface curvature d 2 U/dx 2 (see Eq. 
([15])). Positive (negative) sign of S XtV means that there 
is a material loss (gain) at the crests (d 2 U/dx 2 < 0) and 
a material gain (loss) at the depressions (d 2 U/dx 2 < 0). 
Thus S x . y > (S x, y < 0) translates into a stabilizing 
(destabilizing) action on the surface roughness. Other- 
wise, a positive (negative) sign of S* y means that there 
is a loss (gain) at the crests (d 2 U/dx 2 < 0) and gain 
(loss) at the depressions (d 2 U/dx 2 < 0) of the preferen- 
tially displaced material. Thus this term acts in desta- 
bilizing the composition in all the cases except for the 
situation with no preferential displacement or deposition 
(both alloy components have similar ballistic properties) 
or with flat surfaces (d 2 U/dx 2 = 0). Therefore surface- 
roughness instabilities can start independently from the 
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TABLE I: . Coefficients of the terms of the right-hand side of Eq. (|18[) . Note that the factor p s /A/j at can be absorbed into 
dimensionless diffusivities which would result into a single coefficient in B and B* of cu — 7A 2 /T. 



Y = (yi-Yj) Y* 

v = «W+(l-cM 

D = (Da — Db) Ps/jat D* 

S = [c&D A + (1 - <§)Db] p° s A 2 ^/k B Tj &t B* 



= (l- c %)Yl+c%Y° 

= c° s (l-c? s )v (YZ-Y°) 

= [(1 - c° s )D A + c° s D B ] p%/M 

= c s (l-c s )(D A -D B )p° s A 2 1 /k B Tj &t 



SZ V = S£ v (0) [c° S Y2 + (1 - c&)Yjj] + S^) [4^ eI + (i - 4)s T B el ] 

= co5^ cp + (1 - c°)S d B cp 

S*k n = 4(l-4)\s^ y (6)(Y°~Y^+Sl%(6)(ST-S^)] 

s;f p = co(i - 4)si op - 4(i - Co )^ cp = Co (i - 4)(sl cp - s d /n + ( Co - 4)s£* 

o* c*,ion , o*,dcp /r» 



surface composition effects while surface composition in- 
stabilities cannot occur without coupling to the surface 
roughness. In conclusion, instabilities can occur via sur- 
face roughness effects or via composition-roughness feed- 
back interactions. 

The curvature-dependent deposition terms S%*y and 
S*'y Cp are proportional to the incoming flux j at (Eqs. 
© - (Unj)) • In case that the deposition has a desta- 
bilizing effect, jat is a control parameter. The terms 
~ RS XtV = RS™ + S x *v and ~ R$x,y = -^x!™ + ^x,y P 
explicitly have linear dependencies on R which controls 
the degree of ion-induced effects compared to deposition 
effects. While the exact dependence is determined by 
the details of collision cascades^ 8 - we should note that 
both S x p y (0) and S x c ^ y (9) can change their sign result- 
ing in stabilizing or destabilizing behavior. This can be 
controlled by varying the ion-incidence angled In addi- 
tion, the sign change of S x el (9) and S x p (9) is governed 
by different properties of ion induced collisions cascades 
and thus occur in different angular ranges. 48 As it will 
be shown in the next sections, such general properties of 
Eq. (fT5|) and the dependence of ion-induced effects on R 
and 9 create a rich variety of pattern formations besides 
those induced by deposition. 

As the curvature-dependent terms are fixed, the ratio 
of ion to deposition beam R acts as a control parame- 
ter. If the ion-induced curvature-dependent ballistic term 
combinations (R ^ terms) become negative for a cer- 
tain ion-incidence angle 9 range, their degree of influence 
on the surface dynamics is 'regulated' by R: the higher 
R, the larger the influence of ion effects. If deposition dy- 
namics leads to surface stabilization, then at some critical 
i? cr it instability is expected to result in a growth of a pat- 
tern. In turn, if the deposition dynamics leads to surface 
destabilization, one can find a certain ion-incidence angle 
9 range, where the ion-induced curvature-dependent bal- 
listic term combinations become positive. After reaching 
a critical ion-to-atom arrival ratio i? C rit the growing sur- 
face becomes stabilized. 



B. Stability of modes 

Let us perform a linear stability analysis when the 
pattern growth occurs. For this we consider only a ID 
case along x. The solution of Eq. (fT"8f in the form of 
(U, <fi) = (U°, <fi°) x exp (ikx + rt) represents the growth 
rate r as eigenvalues of the matrix 



A = 



-RS x k 2 - Bk 4 -RY - Dk 2 

-RSlk 2 - B*k 4 -1-RY* -D*k 2 



(19) 



with the eigenvalues 

r = 1 (-/(g) ± Vf(q) 2 -Mq)) , 
f{q) = -trA = B (d + eq + q 2 ) 
a(q) = dot A = / dis q (c + b q + q 2 ) , 

k 2 = q > and the auxiliary quantities 



where 



(20) 



/dis 



BD* - B*D = ^1d a D b > 
B + R(S X D* - S*D + BY* - B*Y) 

/dis 

R [S x + R(S X Y* — YS*)] 



/< 



(.lis 



1 + RY* 
B 

D* + RS X 
B 



> 



(21) 



For a given wave number k, a positive (negative) r value 
indicates that this mode is unstable (stable) and will grow 
(be suppressed) in amplitude. The growth rate r is neg- 
ative (stable) if a(q) > which is fulfilled in the absence 
of deposition and ion irradiation since always /di ss > 0. 
The combination of parameters in /di ss acts as a 'dissi- 
pative' force bringing the system to and maintaining at 
the equilibrium when external influences are absent. In 
order to determine the band of unstable vectors k, we 
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a(q) 



(b) 



f(q) 



i 


a(q)~q(q 2 +bq+c) 


stable 




c>0, b>0 






./''b>-2\lc" / b<-2\|c 




/ / ► 




\J / q 




— ' c>0 b<-2\|c 


Type II 


Type I 



(c) 



f(q)=l+RY + (D +RSx) q + Bq 

stable 
D+RSx>Oy 
Sx>0 > 

Sx<0 




• 0>D+RSx>-21B(l+RY) 



Sx<0 

D+RSx<-2ii(l+RYy 



FIG. 2: (a) The different regimes of unstable modes in terms 
of (b) the parameter a(q) and (c) the parameter f(q) corre- 
sponding to Eq. (|2ip. The arrows indicate possible reachable 
instability regions with continuous change of R. 



proceed first with the stability analysis^ illustrated in 
Fig. [2] . The system is stable for a(q) > and f(q) > 
since then r < for all k and has oscillatory solutions 
if fil) 2 < 4a(g). In the absence of ion irradiation and 
deposition, the condition f(q)>0 is always satisfied (see 
also the discussion above). Within the continuous change 
of the control parameters we can reach the three adjacent 
instability regions from the stable one by the three arrows 
indicated in Fig. [2] (a). The left lower quarter in Fig. [2j 
corresponding to f(q) < and a(q) < 0, can only be 
reached via a special point a(q) = and f(q) — 0. This 



condition requires both a(q) and f(q) to become zero for 
the same value of R and/or deposition S dcp . This can be 
reached only with a unique combination of parameters 
and maximal 2 special wave numbers. This regime does 
not provide any control means and therefore is of low 
relevance to experiments. We are left with two different 
possible paths from stability to instability: (i) a(q) < 
and f(q) > for stationary- growing patterns and (ii) 
f(q) < and a(q) > (f(q)/2) 2 > for oscillatory pat- 
terns since the square-root term becomes purely imagi- 
nary for the growth rate r in Eq. (|2H)) . 



C. Stationary patterns 

Stationary instabilities occur when f(q) stays positive 
and a(q) changes its sign to become negative (see Fig. [2] 
(a)) . This occurs in the range of the two q values 



4 



4c 



(22) 



determined by the sign of parameters b and c which is 
illustrated in Fig. [5] (b) . We see that two real q occur for 
either c < or c> 0, b < — 2x/c. 

Before continuing to discuss the different instability 
types of these two cases lets discuss the physical mean- 
ing of the terms b and c which are crucial for the type 
of instability. The term b has three contributions (Eq. 

my. 

• An R- independent positive term B/f^s ~ jat x 
7 [c° s D A + (1 - c° s )D B ] > related to the film 
growth rate and capillary forces. This term stabi- 
lizes the surface composition and roughness. It in- 
creases concomitantly with the deposition flux j at . 
This term always renders b > for low values of R; 

• A term R(BY* — B*Y) describing how preferential 
sputtering interacts with the curvature to destabi- 
lize the surface. This term is always positive as 
Y* > 0, Y* > Y, B > 0, B > B*. This reflects 
the fact that both surface capillary forces B and to- 
tal sputtering rate Y* act in stabilizing the surface 
roughness and composition, see Eq. ([i"8j); 

• A term R(S X D* — S*D) describing how the 
curvature-dependent (total or preferential) dis- 
placement/deposition processes interact with dif- 
fusion to (de)stabilize the surface. Only this term 
can lead to b < 0, thus to pattern formation. 



The term c in Eq. (f2ij) has two contributions: 
c - R [S x + R{S X Y* - YS*)}. The first one is RS X = 
RS™ + S dcp . This is the total curvature-dependent dis- 
placement/deposition coefficient. The second one repre- 
sents (de) stabilization related to the interaction between 
curvature-dependent displacement/deposition and sput- 
tering (both total and preferential). At R sa (weak 
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or absent ion irradiation), the sign of the curvature- 
dependent deposition term S^ cp determines the slope of 
a(q) growing for small wave numbers q towards negative 
or positive values, and thus determines the instability. 



D. Oscillatory patterns 

When ion irradiation or /and deposition is present, os- 
cillatory behavior might appear if f(q) < or e < —2y/d. 
Such condition yields (see Eqs. (|2"T10 



RS X < - hy/B^ + RY*) + D i 



(23) 



As R > and the expression in the brackets on 
the right side of Eq. (|2U|) is positive, this necessarily 
means S x < and the total curvature-dependent dis- 
placement/deposition term must become destabilizing. 
In addition, this also means that the instability can oc- 
cur only when the control parameter R > exceeds a 
certain critical value. Following the discussion above the 
sign change of f(q) is driven by the interplay between 
the ion-irradiation/deposition induced surface instabil- 
ity and the surface diffusion and preferential sputtering 
rate. This is because the curvature-dependent displace- 
ment/deposition S x < destabilizes while the surface 
diffusion D* > and preferential sputtering Y* > sta- 
bilize the system. 



IV. TYPES OF INSTABILITIES 

In order to understand the interplay of ion-irradiation 
and deposition effects to induce instabilities, we will con- 
centrate separately on the effects produced by ion irra- 
diation and deposition. Their combined effect will be 
discussed at the end. Therefore let us consider first a 
planar growth of the surface (S dep = 0) and concentrate 
exclusively on the instabilities induced by curvature- 
dependent ion displacements. 



A. Ion-irradiation induced instabilities 

1. Instabilities driven by linear terms in q (type II) 

a. Type Ha: curvature- dependent displacement- 
driven pattern formation. First we discuss the case 
when c < 0. This condition can be satisfied for two 
possibilities S x ^ according to Eq. (|2ip . As seen in 
Fig. [5](b), in both cases the instability starts from q = 
and the width of the band Aq of unstable modes is equal 
to the positive solution q£ of Eq. (|22p . For the case of 
S x < 0, expanding in powers of R we get from Eq. (|21[) 
and (1221) the bandwidth of unstable wave vectors 



Afc 2 = Aq 



-R 



S x 
~B 



Near the threshold where a(q) changes the sign it 
becomes small so that r w —a(q)/f(q). The maxi- 
mal growth rate is approximately given by the condition 
da(q)/dq = which leads to 



ftnax = \W b2 - 3c- 6) « ~ R ^§ 



-o{R?). (25) 



The most unstable wavelength diverges when R — > 
while the band width shrinks as ~ VR. This is the so- 
called type-II instability^ 2 . This instability is character- 
istic for ripple formation during ion erosion*^^ Both 
J; max and Afc depend only on the parameters related to 
the surface roughness. Therefore, in this regime (c < 0), 
the pattern formation is driven by ion-induced surface- 
roughening processes. A surface roughness pattern de- 
velops first. The segregation follows due to the h(x,y)- 
c s {x 7 y) coupling. 

Note, that the condition S x < requires some further 
restrictions to render c < 0. If S X Y* — YS* < the 
instability condition c < is fulfilled for all R . If S X Y* — 
YS* > the ion-to-atom arrival ratio should not exceed 
a critical value R < R C nt.i with 



R 



crit.l 



S X Y* - S*Y 



(26) 



o(R 2 ). 



(24) 



b. Type lib: preferential curvature-dependent dis- 
placement and sputtering. For the case c < 0, let us con- 
sider now the possibility S x > 0, i.e. the total curvature- 
dependent displacement acts in stabilizing the surface. 
Such conditions requires S X Y* — YS* < (see accord- 
ing to Eq. (pT]l. As both S x and Y* (total sputtering 
coefficient) are positive, the instability is driven by the 
combination YS* > 0, i.e. by preferential sputtering and 
preferential curvature-dependent displacement. The con- 
dition c < requires that these both terms must have the 
same sign (S* > as Y* > by definition). 

This represents the following scenario: 

• the element segregation occurs on the surface as the 
preferentially displaced element is depleted from 
the crests and enriched in the depressions (term 

inEq. (UK, 

• preferential sputtering of the same element occurs 
(term -RY<j> in Eq. ([j"8]t. 

• as this element tends to segregate at the depres- 
sions, preferential sputtering of the depressions 
takes place which results in the increase of the sur- 
face roughness, i.e. an instability occurs. 

The instability is also of type II following the classifi- 
cation by Cross&Greenside^. The main difference to the 
type II a instability is that the ion-to-atom arrival ratio 
R must be larger than a critical value Eq. (|2l)]) . Close to 
the instability threshold one can write R = i? C rit,i + S. 
Introducing this expansion into Eq. (|22[) and leaving only 
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linear terms in 6 to get Ac = —S x /fdis8 and the following 
expression for the band width Aq 

Ak 2 = Aq « -j^-S (27) 

/disherit 

where 6 cr it is the value of the parameter b at i? C rit,i- The 
bandwidth shrinks again as Afc ~ y/S. The maximal 
growth appears again as in Eq. (|24j) and Eq. ([25]) nearly 
at the minimum of a(q) leading to fc^ ax = g max ~ -g 2 as 
above. 

Note, that the destabilizing effect is non-linear R, i.e. 
~ R 2 . This is due to the fact that both the instability- 
inducing terms originate from ion irradiation. Both 
terms are based on the differences in ballistic properties 
of the alloy constituents. However, one term refers to 
general alloy sputtering properties (preferential sputter- 
ing) while the other refers to a topographic effect (pref- 
erential curvature-dependent displacement). Note that 
their product comes from multiplication of off-diagonal 
terms of A in Eq. (fTO)) . Off-diagonal elements describe 
the feedback interactions between the surface rough- 
ness and the composition. This brings us to the two 
conclusions: (i) instabilities can be generated by ion- 
induced h(x,y)-cs(x,y) feedback interactions and (ii) the 
strength of the ion-induced feedback interactions is driven 
by differences of material ballistic properties. The larger 
the difference in ballistic properties (sputter/relocation 
yields and their sensitivity to the surface curvature), the 
larger is the probability for this effect to occur. 

2. Instabilities driven by quadratic terms in q (type I) 

The term a(q) can become negative not only due to lin- 
ear terms in q proportional to c but also due to quadratic 
terms in q proportional to b (see Eq. (j2"0")) '). Let us con- 
sider the case when c > and b < —2^/c. The slope 
of a(q) is positive close to zero due to the positive lin- 
ear term cq. Only for larger q values the term ~ — \b\q 2 
counteracts the positive increase leading to a(q) < (see 
middle plot (b) of Fig. [5]). This results in a band of 
unstable q vectors appearing in a narrow q range with 
q > 0. Therefore, close to the instability threshold the 
wavelength is not diverging. The instability will not oc- 
cur for any small R, but R needs to reach a critical value 

#crit,2 given by b 2 = 4c (28) 

from Eq. (|2"Tj) . 

The condition b < requires S*D > (sec Eq. (j2Tf ). 
The terms S* and D describe the preferential curvature- 
dependent displacement rate and preferential diffusion 
of the alloy constituents, respectively (see Table HJ. This 
represents the following scenario: 

• for positive (negative) S* values the preferentially 
displaced element is depleted from the crests (de- 
pressions) and accumulates in the depressions (on 
the crests) (the term RS*^- in Eq. ([15]). 



• such a segregation induces a concentration gradient 
between crests and depressions, 

• a diffusion-driven transport sets in to homogenize 
the element distribution on the surface (the term 
DV 2 <j) in Eq. (p]). 

• the element with a larger diffusion coefficient dif- 
fuses from the depressions to the crests, 

• as this element tends to segregate at depressions, 
preferential diffusion of this element from depres- 
sions results in the increase of the surface rough- 
ness, i.e. an instability occurs. 

Let us assume that close to the instability threshold 
R = -R C rit-2+<5 and therefore (6/2) 2 — c = const X S+... > 
0. According to Eq. (|2"2"1) . the band width takes the form 
Aq = % — % = \/b 2 — 4c yielding 

Ak x = yfAq~~8 1/4 . (29) 

The band center can be determined from the condition 
da(q)/dq = 0. Leaving only linear terms close to the 
instability results in 



kmax » \ --bait + 7T f X 5. (30) 

V O I & crit | 

The smoothing effect of the curvature-dependent dis- 
placement coefficient S and diffusion D* separately sta- 
bilize the surface and the composition at small and large 
wave numbers, respectively. This allows to obtain a nar- 
row band of unstable wave numbers, i.e. the formation of 
an ordered structure. According to Cross&Greenside^ 2 
this is characteristic for the so-called Type-I instability. 
Similar effect has been predicted for ion-erosion of al- 
loy surfaces.— In contrast to the ion-erosion case where 
the surface atoms are always exposed to ion irradiation, 
here the interplay of the ion-irradiation and diffusivity is 
'damped' by the depositing species (both parameters R 
and D are divided by j a t , see Table Q} . 

Note, that the destabilizing effect comes as a linear ef- 
fect in R. This is due to the fact that only one term in the 
product S*D refers to ion- induced effects (preferential 
curvature-dependent displacement). Similarly to the in- 
stability type lib, the product of terms S*D inducing this 
type of instability originates from multiplication of off- 
diagonal terms of A in Eq. (|19[) which describe h(x,y)- 
c s{x,y) feedback interactions. The preferentially dis- 
placed element must also have a larger diffusivity to keep 
S*D > 0. The above discussion leads to the following 
conclusion: external preferential ion-induced curvature- 
dependent processes and material-inherent diffusive pro- 
cesses can couple to induce an instability. The larger 
the difference in curvature-dependent ballistic and sur- 
face diffusion properties, the larger is the probability for 
this type of instability to occur. 
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3. Special paths from type-II to type-I instabilities 

By increasing R one can observe a change from type-II 
instability to a stable behavior and further to the ap- 
pearance of type-I instability. Such a behavior requires 
that the total curvature-dependent displacement coef- 
ficient has a destabilizing effect on the surface rough- 
ness (S x < 0) and a stabilizing effect due to synergistic 
effects of the preferential curvature-dependent displace- 
ment rate and the preferential sputtering, S X Y* — S*Y > 
0. 

Mathematically, this is due to the fact that by in- 
creasing the ion-to-atom arrival ratio R, the coefficient 
c - R[S X + R(S X Y* - YS*)} (first order in q term) can 
change the sign from negative (first term of c) to positive 
(second term of c). Hence we have the transition: type- 
II instability— > stable behavior. With further increase of 
R, the coefficient b > (second order in q term) might 
change the sign, which results in type-I instability. 

The further discussion of these parameters is quite 
complex and is outside of the scope of the present paper. 
Our purpose is only to point out here that certain alloy 
deposition conditions can lead to type-II— ^stable— >type-I 
transitions when increasing the ion-to-atom arrival ratio 
R. 

Jj.. Oscillatory instability 

The instability growth exponent r can also change the 
sign when /(g) changes the sign from positive to negative 
while a(q) remains positive (see Fig. [5] (a)). This leads 
to an increase of the instability amplitude due to the 
first term r ~ ~f(o) > 0. However, the second term of 
r ~ ■J ' f{q) 2 — 4a(q) (Eq. (|20|) remains imaginary, hence 
we have an oscillating behavior with the increasing am- 
plitude. The necessary condition follows from the expres- 
sion of f(q) in Eq. 

e < -2Vd. (31) 

For e < 0, the total curvature-dependent displacement 
yields S x must be negative, thus destabilizing (see Eqs. 
([20]) and ([2~Tj) ). Like the type-II a instability, it is also 
driven by ion-induced surface-roughness processes. 

The oscillatory instability Re(r) > requires exactly 
the same g-range where /(g) becomes negative. The in- 
stability occurs when R surpasses a critical value 

R CI it, 3 defined by e 2 = 4d. (32) 

Near the critical threshold i? C rit,3 we can expand again 
e 2 /4— d rs const x 6 and obtain the unstable wave number 
band width in lowest order 

Afc 2 = Aq-5 1 / 2 . (33) 

The maximal growth rate is given now approximately 
by the maximum of /(g) which leads to 

W — <Zmax — 7J ~ 2 ~2B 



with the value e crit = (D* + R crit)3 S x ) /B (Eq. ([2l|)l. 
Note that the characteristics of this type of instability 
differs from those of the oscillatory type-Ill instability 
from Cross&Greenside£2. The real part of r is not zero, 
therefore the maximum growth rate occurs at a finite 
wavelength in Eq. (|34j) . At the onset of the instability 
where R = i? C rit,3 we have Re (r) = —f(q) = 0, and the 
instability is purely oscillatory. Therefore the type-Ill of 
Cross&Greenside is a special case of our type I in the 
sense of finite q values as an onset. 

In more detail, additionally to the condition e < 
—2yd, two more conditions must be satisfied to render 
a(q) > 0, namely 

c> 0; b 2 < 4c. (35) 

The condition S x < means that c < for small R which 
contradicts the necessary conditions Eq. ([33)) for the os- 
cillating instability to occur. If the parameter c would 
not change its sign to positive with increasing R, then 
the change of sign of f(q) would result only in the change 
of the instability type, not in the change from stable to 
unstable situation. We are limiting our discussion only to 
the transition stable— > unstable. Therefore we consider 
that c < ^ c> by increasing R. In addition, the 
change in the sign of c must occur for R < R cr it,3- 

On one hand, after reaching this stabilizing condition, 
a further increase in R keeps a(q) > since both lin- 
ear and quadratic terms proportional to c and b, respec- 
tively, are positive. On the other hand, it induces the 
sign change for /(g), thus the occurrence of the insta- 
bility. Recalling Eq. (|2"3"|) . one can see that the stabi- 
lizing effects are all proportional to diffusion (both D 
and ~ \f~B ~ VD). Therefore, for the transition type- 
IIa — >stable— ^oscillating instability to occur, the diffu- 
sivity must be large enough to allow stabilization of the 
surface while increasing the ion-to-atom arrival ratio R 
(transition II a — ^stable). Then further increase in R un- 
avoidably leads to the change of sign of f(q) and hence 
the transition stable— >oscillating instability. Physically, 
such a control over surface diffusivity can be achieved 
independently from the ion-induced effects by changing 
the substrate temperature. 

The imaginary part of the instability growth rate r is 
given by Im(r) = ^/4a(g) — f(q) 2 - It defines the oscilla- 
tion frequency in thickness units. For the growing film 
this means that for a given lateral position (x, y), the con- 
centration will oscillate with time. Since the film is grow- 
ing, this translates into an oscillatory concentration with 
thickness. Not only a laterally ordered nanocomposite 
structure forms, but also vertical composition modula- 
tions occur with the period ps 1 /Im(r) . Such an oscillat- 
ing instability establishes a self-organized 3D-multilayer 
structure. This demonstrates that ion-irradiated surfaces 
during alloy film growth can induce not only lateral but 
also vertical periodic structures. 
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B. Deposition induced instabilities 

Now we investigate the possibilities to reach the in- 
stable regions in Fig. [5] due to dynamic processes in- 
duced by deposition only. In the absence of ion irra- 
diation, the terms c and b become c = S p //di s and 
b = (B + S dcp D* - S* dcp D)/f dis (Eq. (EI])). Stationary 
patterns occur for wave numbers in the interval Eq. ([22]) 
and for c < or c > 0, b < —lyfc. This translates into 
two types of instability 

Typella :c < O 5 dcp < 
Typel :c> 0,6 < -2«/c 

O 5 dcp > 0, 5 dcp L>* - S* dep D < —B — 2f dis v f S d ° p . 

(36) 

This shows that for reaching the type-I instability a suf- 
ficient preferential deposition S^ p — < S c is neces- 
sary. We see that for reaching the unstable region of the 
lower right quarter in Fig. [5^ only type I and type Ha 
instabilities are possible from the ones identified above. 

Recall, that the the upper left quarter of Fig. [5] (a) 
represents the oscillatory instability of type III above. 
For such oscillatory patterns we have e = (D* + S dep )/B 
and <i=l/£>>0in the absence of ion irradiation R = 0. 
f(q) has to change the sign and a(q) > which translates 
into the conditions of Eq. ([33)) . The condition c > 
becomes S dcp > and the second condition of Eq. ([33)) 
becomes D* + S dcp < —2\fB. This is excluded since 
the left side is positive. We conclude that no oscillatory 
instability can be reached at perpendicular atom incidence 
due to curvature- dependent deposition effects. Similarly, 
it was found by Guyer and Voorhees (Fig. 4 in Ref. I33T ) 
that the imaginary part of the instability growth rate is 
connected with a negative, i.e. damping, real part. For 
this type of instability to occur one needs oblique angle 
deposition^ 

Since in the absence of ion irradiation (R = 0) no oscil- 
lating instability can be induced due to the deposition ef- 
fects we can limit ourselves only to the deposition induced 
type-I and type-II instabilities. These are controlled by 
the terms linear and quadratic in q proportional to c and 
6, respectively, as discussed above. 

Similarly to the ion-induced instability, the sum of the 
dynamic deposition effects must be destabilizing to pro- 
duce the type-II Q instability. This yields that crests are 
growing faster than depressions. The band width and the 
fastest growing wave number are (Eqs. (IM)) and [25]) 

cdep 



s dcp 

^max ~ ^g - ~ Jat- (37) 

The fastest growing wavelength is diverging for small 
fluxes. The instability growth is the interplay between 
the deposition dynamics and the diffusive surface relax- 
ation. The strength of deposition effects increases con- 



comitantly with the deposition flux. The surface dif- 
fusivity is constant for a given substrate temperature. 
The maximum wavelength (wave number) decreases (in- 
creases) with the deposition flux since for higher deposi- 
tion rate the instability-causing fluxes can compete with 
diffusive-relaxation fluxes at smaller scales. 

Again, similarly to curvature-dependent ion-induced 
instabilities, the type-I instability requires (i) a suffi- 
cient level of curvature-dependent preferential deposi- 
tion — < S c and (ii) a coupling between the 
curvature-dependent preferential deposition and the pref- 
erential diffusivity. The curvature-dependent preferential 
deposition means that one element accumulates prefer- 
entially at the crests or at the depressions. This is a 
result of different dynamic interaction strengths of alloy 
components with the topographic features and defects or 
differences in latent mobilities. The element with larger 
latent mobility is expected to funnel down more easily. 
The element with stronger dynamic interactions with the 
surfaces is expected to be 'steered' towards the crests. 
The following scenario occurs: 

• for positive (negative) S* values, the element with 
stronger (weaker) interaction strength at the sur- 
face is depleted from the crests (depressions) and 
accumulated in the depressions (on the crests), 

• such a segregation induces a concentration gradient 
between crests and depressions, 

• a diffusion-driven transport sets in to homogenize 
the element distribution on the surface, 

• the element with a larger diffusion coefficient dif- 
fuses from the depressions to the crests, 

• again, as this element tends to segregate at depres- 
sions, preferential diffusion of this element from 
depressions results in the increase of the surface 
roughness, i.e. an instability occurs. 

For this effect to occur, a critical deposition flux j a t,crit 
must be reached. Near the instability threshold one can 
assume j at = J a t,crit + e an d the band width takes the 
form 

Ak x ~ e 1/4 . (38) 

The fastest growing wavelength has an identical expres- 
sion as the one obtained for ion-induced instabilities (Eq. 
(|30[) with 6 cr j t defined at j a t,crit)- Similarly to the ion- 
induced instability case, the smoothing effect of the de- 
position dynamics and diffusion separately stabilize the 
surface and the composition at small and large wave num- 
bers, respectively. This allows to obtain a narrow band of 
unstable wave numbers, i.e. the formation of an ordered 
structure. 
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C. Interplay of ion-irradiation and deposition 
induced instabilities 

Curvature-dependent ion-induced displacements and 
deposition processes can act in the same direction (sta- 
bilizing or destabilizing) or in opposite directions (one is 
stabilizing, another destabilizing). If both processes are 
destabilizing, then all the conclusions made in the sec- 
tion |IVX] are valid. Mathematically this means that the 
expression for the curvature-dependent instability driv- 
ing coefficients S x and S* must contain terms related to 
ion-induced and deposition-induced effects. Physically it 
means that instabilities are controlled by both the atomic 
flux j at and the ion-to-atom ratio R. The critical values 
to induce type-IIf,, I and III instabilities will also depend 
on both of these parameters. 

In the case that the curvature-dependent ion-induced 
and deposition-induced processes act in opposite direc- 
tions, the situation differs. For a given atomic flux, a 
critical i? cr ;t will exist for each type of instability, also 
for type-II a . This is due to the fact that ions (deposit- 
ing atoms) have to compensate not only the relaxation 
effects induced by diffusive fluxes, but also induced by de- 
positing atoms (ions). For the type-II a instability, which 
occurs at any small R or j at values in the cases consid- 
ered above, the critical value is defined by the condition 
c = 0. For other types of instabilities the mathematical 
expressions presented in Section TlV Al are valid again, but 
one has always to consider that critical ion-to-atom ratio 
values must be larger. 

Note that the stabilizing or destabilizing character of 
ion irradiation is controlled by the ion incidence angle 
9. This angle can be adjusted independently from the 
deposition fluxes. If deposition is destabilizing the film 
surface, and a homogeneous smooth film growth is de- 
sired, the ion incidence angle can be adjusted to produce 
a stabilizing effect. Then above the critical ion flux value, 
the ion effects compensate deposition-induced instabili- 
ties. In the opposite case when a heterogeneous structure 
is desired, the ion-incidence angle 9 must be adjusted to 
act in a destabilizing manner. Though the diffusive and 
deposition-induced fluxes act in a stabilizing manner, a 
critical ion flux can be reached when an instability oc- 
curs. Similarly, if ion irradiation is used to assist the 
growing film and its effect is destabilizing, an increase 
in stabilizing the atomic flux j at can result in a decrease 
of the ion-to-atom ratio R below the critical value. Ho- 
mogeneous and smooth growth is produced. This leads 
to the following conclusion: ion irradiation (deposition) 
can be used to stabilize the surface roughness and com- 
position if destabilization occurs due to deposition (ion 
irradiation) effects. 

The surface curvature affects not only instabilities, but 
also the deposition/re-sputtering balance. This balance 
is controlled by the terms proportional to sputtering and 
deposition coefficients in the total curvature-dependent 
rate Sx (see Table HI . If dynamic deposition effects result 
in matter accommodation on (depletion from) some to- 



pographic features where sputtering is larger (smaller), 
a decrease (increase) of the total net growth rate V de- 
fined by Eq. (0} results. In a similar manner, surface 
curvature can affect the steady-state surface concentra- 
tion Cg. The balance of the species is regulated by the 
deposition/sputtering terms in the preferential rate S*. 

If a certain element is preferentially deposited (de- 
pleted) on certain topographic features where its sput- 
tering rate is larger (smaller), a net decrease (increase) 
of this element happens in comparison to the steady-state 
value c° s defined by Eq. ([5]) . Also any local enrichment of 
an clement results in the total loss due to the sputtering 
term —Y*<f> (see Eq. (jTSj)) . Note that this balance can be 
regulated by the ion incidence angle 9 which can change 
and even inverse the behavior of sputtering for different 
topographic features^,. This brings us to the conclusion 
that surface topography does not only influence the local 
redistribution of the deposited matter but also influences 
the net growth rate and the average steady state surface 
composition. 

V. PHASE DIAGRAM OF EXPERIMENTAL 
PARAMETERS 

A. Summary on boundary conditions 

In order to discuss the phase diagram of instabilities 
induced by curvature-dependent deposition and ion irra- 
diation effects we proceed in two steps. First, we discuss 
the effects of ion irradiation, i.e. S dcp — exclusively. 
This will provide us with a unique set of three parame- 
ters which will be expressed as combinations of sputtering 
and relocation yields. Such parameters will give some in- 
sights in the internal structure of the instability driving 
curvature-dependent terms S x °y and S*)° n . In a second 
step we introduce the deposition terms. 

1. Invariant measures 

Now we will link all conditions discussed so far to the 
original parameters in Table U in a transparent way. To 
this aim we use three positive combinations of relocation 
and sputtering, 
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and rewrite the conditions by equivalent ones as follows 
S ion Y * _ y S *ion ^ ^ s s x p (9)§-Sf(9)P 3 , 

S™ ^ o S S J>(6) ^ Sf(9)P 2 . (40) 
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One can understand Pi (P 2 ) as the ratio of preferential 
(total) relocation to preferential (total) sputtering. P 3 is 
the sum of species-dependent relocation and sputtering 
coefficient ratios. The parameters Pi,P3 are controlling 
parameters of the surface concentration while the param- 
eter P2 is related to the local height. 

Interestingly, these parameters are bounded by 

Pi ^ P 2 ^ P 3 o Y°Sf ^ YsSt (41) 

and no other possibilities. Here we have assumed that 
the preferential sputtered species A is also preferentially 
relocated, so that P\ > 0. 

The decisive condition Y%S% 1 ^ YgS 1 ^ in Eq. flU]) 
may seem to be only determined by ballistic properties 
of given materials A and B, i.e. ratios of sputtering and 
relocation yields. Therefore low level of control from the 
experimental parameters (ion energy, type, etc..) may 
be expected. However, one should note that varying the 
ion energy close to the sputtering threshold, the relo- 
cation yield varies slowly in contrary to the sputtering 
yield. The latter is zero below the threshold, indepen- 
dently of the relocation yield, and rises to some finite 
value above it. This allows to influence significantly the 
ratios Y® / to match the corresponding requirements 
for instabilities. 

With Eq. (|40|) the conditions for the instabilities can 
now be translated into relations between sputtering and 
relocation. 



2. Ranges of instability types 

a. Type I This type of instability requires b < as 
the first necessary condition. The exact one, b < — 
is ensured by exceeding the critical value R > P C rit,2 of 
Eq. (|2"5)) . From the definition of b in Eq. (|2"T|) we see that 
this means S X D* - S*D < B*Y - BY* < 0. The last 
inequality follows from Table U which gives BY* — B*Y = 
[D B (1 - c° s )Y A + c° s D A Y B }p s A 2 1 /k B Tj at > 0. Further- 
more, B* ^ always if D ^ which means preferential 
diffusion D A D B . Consequently, the condition b < 
translates into 

S x < 0,S*^0 for D A ^D B . (42) 

The second necessary condition c > with S x < re- 
quires S X Y* - YS* > and R > P cr it,i- Using Eq. (gDJ) 
this translates into 

Sf(9) ^ : 

-Sf(9)P 2 >S x p >-Sf(9)P 3 for H^ll 
Sf(9)^0,D A >D B : 

-S x c \9)P 2 >S sp >-S x c \9)P 1 for Hr^||. 

(43) 



It leads us to the conclusion that S s x p (9) and S r x e \9) 
must have opposite signs. In the case sputtering acts 
in stabilizing the surface roughness S x p (9) > 0, the re- 
locations must act to destabilize it (0) < 0, or the 
other way around. Therefore the instability can occur 
only within the ion incidence angle range where sputter- 
ing and relocation have opposite impact on the system 
stability. As mentioned before, the signs of S x cl (9) and 
S x p (9) is governed by different properties of ion-induced 
collisions cascades and thus occur in different angular 
ranges 42, Therefore there are large angular ranges where 
these two quantities have opposite signs 4^ 

Note that relocation yields are usually much larger 
than sputtering yields, Sf 1 >> S sp , thus Pi,P 3 >> 1. 
Maximum values of S x p (9) and S%(9) are of the order of 
unity4£ This means that not only the signs of S x cl (9) and 
S x p (9) must be opposite but also |S* el (0)| << \S sp (9)\. 
The condition of the relocation term being close to zero 
confines the ion incidence angles around ^-values where 
the relocation term changes its sign, thus its surface sta- 
bilization behavior. This brings to the conclusion that 
the type-I instability is confined to small angular ranges 
where S x cl {9) w 0. 

b. Type Ha From the discussion above follows, 
that the curvature-dependent displacement-driven pat- 
tern formation is possible for S x < and S X Y*—YS* < 
for all R or S X Y* -YS* > but R < R aitA via Eq. (JH]). 
Therefore we have according to Eq. (l40l) the following 
possibilities 

Sf(9)>0,VR:S x p <-Sf(9)P 2 for ^ ^ || 

^A a B 

Sf(9)<0yR:S sp <-Sf(9)P : , for |J $ || 

2 ^A S B 

^ el (^) §S 0, P < P C rit,l : 

- Sf(9)P 3 < S sp < -Sf(9)P 2 for H^ll 

S A S B 

(44) 

where we abbreviate the equations such that the up- 
per/lower index of P corresponds to the upper/lower 
inequality. As expected, this type of instability is not 
related to Pi which contains in its expression the terms 
related to preferential relocation or sputtering ( see Eq. 
©). 

c. Type lib Following Eq. (|40l) we have as condi- 
tion for the preferential curvature-dependent displace- 
ment and sputtering, S x > and S X Y* - YS* < 0. This 
leads to 

Sf(6) ^0,R>R crit , 1 : 

-Sf(9)P 2 <S sp <-Sf(9)P 3 for Ui^H. 

a A d b 

(45) 

Here we abbreviate again the equations such that the 
upper/lower inequalities correspond to each other. Al- 
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though the preferential sputtering is necessary for this 
type of instability (see Section IIV Al) , it does not appear 
explicitly. The instability is governed by differences in 
ballistic property ratios of each component. 

d. Type III The necessary conditions for this type 
instability are e < which means that S x < and c > 0. 
With the analogous discussion as above, this translates 
into the condition 



Sf(6) ^ : 

- Sf{9)P 2 > S*> > Sf(9)P 3 for 
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(46) 



B. Phase diagram of instabilities due to 
curvature-dependent ion irradiation effects 



We collect the different identified areas of Eqs. (14"31) . 
(gU), (g5]), (gU) together in Fig. [3J The diagram axis are 
S x cl (9) and S x p (9). The ion-incidence angle 9 is noted 
here intentionally to emphasize that it controls the sign 
of the curvature-dependent relocation and sputtering co- 
efficients. Note again that S x el (9) and S x p (9) change their 
sign at different ion-incidence angles 9.— From the dis- 
cussion above the relative sign changes of sputtering and 
relocation decides which branch of instability type we 
have to expect. Physically, different areas of this dia- 
gram could be reached by varying the incident ion angle, 
the energy or the species type. Three critical values of 
i? cr it discussed above will determine where the relevant 
instability types can occur. 
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Let us first consider the case < -^t- 
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We see that 

the major range is given by the type Ha instability which 
is all the area below the —S^Pz line. Above the — S x cl Pi 
line for S x cl > and above the ~S x cl P 3 line for S 1 / < 
this instability occurs only for a restrictive range of R as 
indicated in the figure. The type lib instability occurs 
only for S x cl > and between the -S r x cl P 2 and -S x el P 3 
lines. The type I instability can occur for both signs of 
S x el . The type III oscillating instability only for negative 
S x cl < between the -S 1 x cl P 2 and -S T x el P 3 lines. The 
overlapping regions in the figure are separated by critical 
values of R via Eqs. flU, (gSJ), ([32]). 



The analogous figure for the case 



> 



Y" 

B 



is plotted 



in Fig. [3J below. One sees that the regions are mirrored 
at the line Si el = Si p . 



C. Phase diagram of instabilities due to combined 
curvature-dependent ion irradiation and deposition 
effects 

Now we include the deposition terms by replacing for- 
mally 
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in the definition of the invariant measures Eq. (135)) . This 
accounts for the deposition-dependent parameters in ta- 
ble HJ We have to change therefore 
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where analogously to Eq. (I4T1) 
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(50) 



The invariant measure Q3 represents the ratio between 
individual curvature-dependent deposition and sputter- 
ing weighted with the incoming flux composition. The 
parameter Q\ is the ratio of curvature-dependent prefer- 
ential deposition to preferential sputtering. The param- 
eter Q 2 represents the ratio of total curvature-dependent 
deposition and sputtering rates. The parameters Qi,Q 3 
are controlling parameters of the surface concentration 
while the parameter Q 2 is related to the local height. 
The two possible orderings of Qi in Eq. (l50l) can be seen 
as preferential curvature-depedent deposition to prefer- 
ential sputtering ratio of species A. 

We see from Eq. (|4"3j) that the matter deposition shifts 
all straight lines in Fig. [3J by a constant Qi/R. There- 
fore, the different cases discussed in Fig. |3J start to over- 
lap and form a new area of phase diagram as seen in 
Fig. 2J This overlap region must contain both condi- 
tions from overlapping areas. Choosing first the case of 
preferential deposition to sputtering ratio of species A, 
> Qi > Q 2 > Q 3 , and the situation of preferential 
relocation to sputtering ratio of species A in Fig. [3J (a) , 
the new overlapping region formed by a triangle contains 
two phases (Fig. 2]) 

type/ with R > max{R CI it,i, RcAtp), D A > D B or 
type Ila with R < R CT i ty i (51) 



16 



and no oscillatory behavior which is restricted to the 
left upper region exclusively due to ion-deposition effects. 
The case of Q\ > Q2 > Q3 > would shift all lines of 
Fig. 2] downwards without any qualitative change of the 
discussion. The relative shift of the axes are controlled 
again by the order parameter R. 

Concerning the case of preferential curvature- 
dependent deposition to sputtering ratio of species B, 
Qi < Qi < Q3, together with the preferential relocation 
to sputtering ratio of species A in Fig. [3] (a), it turns out 
that the new triangle region has to be an overlap of type 
II a and type lit, instability which is impossible. There- 
fore in this case the triangle region does not specify any 
instability. 

The opposite case of preferential relocation to sputter- 
ing ratio of species B in Fig. |3] (b) just mirrors Fig. [3] (a) 
on the axes S* p (9) — S 1 x cl (9) and all discussions above 
apply accordingly. 

The incident angle of ion irradiation determines the ac- 
tual values of the curvature-dependent sputtering S^P(8) 
and relocation S^ cl (9) but not the P, slopes of the lines. 
The angle of deposition determines the actual value of 
the Sf op parameters and therefore the Qi offsets. One 
sees that the actual angle dependencies would determine 
which of the cases apply. The striking observation is that 
despite the undetermined and complicated angular depen- 
dence of these parameters one can give a phase transition 
with all possible occurring instabilities. 



D. Relevance to experimental film growth 
conditions 

The linear analysis above demonstrates that ion- 
irradiation during the alloy film growth is expected to 
play a crucial role on the stability or instability of the 
growing surface. It is controlled by the ion-to-atom ar- 
rival ratio R which can be varied in a broad range suited 
for numerous experimental conditions. The ion incidence 
angle 9 determines the sign of S sp (9) and S rel {9) and 
thus whether and in which way ions act in destabilizing 
(stabilizing) the surface roughness and composition. 

Experiments show that energetic ion assistance can 
act in suppressing the Stranki-Krastanov growth mode of 
SiGe alloy films and stabilizing the smooth growth— An 
ion-induced stabilizing effect has been also demonstrated 
for C-Ti nanocomposite films^— Moreover, it has been 
experimentally varified for carbon films that the assist- 
ing ion irradiation during ion-beam assisted deposition 
(IBAD) can cause the formation of the periodic surface 
roughness patterns—"— These experimental findings are 
in qualitatively agreement with the statement on the im- 
portance of ions to control the growing surface morphol- 
ogy. However, the lack of detailed dependencies of the 
surface nanopattern structure on the ion-to-atom arrival 
ratio R does not allow any quantitative comparison with 
the results of this study. 

Nevertheless, this shows the potential of ion-beam- 



assisted deposition (IBAD) as a powerful tool to investi- 
gate the patterning during the film growth. In contrast 
to ion-erosion case, it allows to produce nanostructured 
surfaces whose material is different from that of the sub- 
strate. In IBAD the deposition is decoupled from the 
assisting ion bear m 40 ' 41 Therefore the control by means 
of R and 9 is possible. The theory developed in this work 
shows that these parameters not only induce the forma- 
tion of a surface pattern, but can also control their struc- 
tural properties. Further experimental work is needed to 
reveal this potential. 

From the theory point of view, it has been demon- 
strated here that ion irradiation together with matter 
deposition can induce phase separation in alloys during 
ion-assisted film growth. Therefore it has similar effect as 
the substrate misfit and compositional stresses^ - — Sim- 
ilarly to the stress effects, the ion-induced compositional 
pattern couples to the surface roughness. Composition- 
ally modulated surface roughness patterns occur. Their 
behavior close to the instability threshold depends on 
the regime: driven by ion-induced surface roughness pro- 
cesses (types II a and III) or by roughness-composition 
feedback interactions (types I and lib). 

Of particular interest is the possibility to produce pat- 
terns with a high degree of ordering. Our theory predicts 
that if the conditions for the type I instability are satis- 
fied, only a narrow band of unstable wave vectors is pro- 
duced. This regime is driven by the coupling of preferen- 
tial curvature-dependent displacement/deposition with 
preferential diffusivity. The conditions require that the 
preferentially displaced material should also be the one 
with a higher diffusion coefficient. This is not always the 
case for thermal diffusion. However, at low deposition 
temperatures ion-irradiation effects on the surface diffu- 
sivity can become substantial ! 71 ' 72 This is expected to act 
on the preferentially displaced element and, therefore, to 
favor the occurrence of this type of instability. 

On one hand, this opens an alternative opportunity 
to produce highly compositionally ordered surfaces. On 
the other hand, such nanopattern is constantly buried 
by randomly depositing species. The formation of a 3D 
nanopattern occurs with an ordered heterogeneous struc- 
ture, hence the control of the resulting nanocomposite 
morphology. The selected wavelength will be finally de- 
termined by non-linear terms stabilizing the initial ex- 
ponential growth of unstable modes. The tilt of such 
structures in relation to the film surface will be propor- 
tional to the drift terms v ■ dh/dx and v* ■ dh/dx (see 
Eq. (|18|l) and thus determined by R. 

In the Si ion-erosion case, it has been recently 
demonstrated that sputtering is not needed to induce 
instabilities! 49 ' 50 While the effect remains to be proven 
for other material systems, this finding shows that one 
can work in the regime below or close to the sputter- 
ing threshold. Sputtering is to be avoided during the 
growth due to the loss of material. It limits the ion-to- 
atom arrival ratio to induce instabilities as complete re- 
sputtering might occur earlier than any instabilities are 
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induced. Besides, for ions exhibiting energies below the 
sputtering threshold, film composition is the same as that 
of the incoming depositing flux j at . Following the theory 
outlined in this work, the ion-induced instabilities are 
only driven by relocation terms. Assuming a simple an- 
gular dependency for the relocation term ~ 008(26*) from 
the literature , 47 ' 48 a relocation driven pattern growth is 
expected for ion incidences 9 >~ 45°. If alloy growth is 
desired without formation of any secondary phases, one 
should grow below the sputtering threshold and at low 
ion-incidence angles 9. 

Note that the above model concerns experimental sit- 
uations where surface displacements prevail over bulk 
displacements^ It is observed experimentally that ion 
irradiation also can induce ballistic transport within a 
'bulk' of a confined thin sub-surface layer. This results in 
vertical compositional patternsi^i 7 ^— For comparable 
fractions of bulk and surface displacements, both mecha- 
nisms - ion induced layering and surface instabilities - are 
expected to act simultaneously, and the produced struc- 
ture will be the result of the competition between both 
effects. Our results indicate that vertical layering can 
also occur only due to surface effects when the oscillat- 
ing instability occurs. 



VI. CONCLUSIONS 

Linear instability analysis has been carried out to 
study the opportunities to induce pattern formation dur- 
ing ion-assisted alloy growth. Curvature-dependent ion- 
induced relocation and sputtering effects have been con- 
sidered in the model as well as slope and curvature- 
dependent dynamic deposition effects. 

It has been demonstrated that during alloy film growth 
energetic ions (i) can induce phase separation and (ii) 
drive the system towards stationary or oscillating sur- 
face pattern formation. As the surface roughness couples 
to the composition, this results in compositionally mod- 
ulated roughness patterns. The control parameter for 
'switching' and 'amplifying' such instabilities is identified 
to be the ion-to-atom arrival ratio R. The ion-incidence 
angle 9 determines whether and in which way the ions 
have a destabilizing effect on the surface roughness and 
composition. 

Three different regimes were found close to the thresh- 
old driven by surface roughness processes or composition- 
roughness feedback interactions. The composition- 
roughness feedback interactions are governed by differ- 
ences in ballistic, dynamic and thermal properties of the 
alloy constituents. Two types of stationary instabili- 
ties and an oscillatory one have been identified which 
correspond to the classification presented by Cross and 
Greenside£2: 

• Type-I instability which is driven by synergistic in- 
teractions of preferential diffusivity and preferential 
curvature-dependent displacement / deposition rate. 



It is characterized by a narrow unstable wavelength 
band starting at finite wavelength whose width Afc 
varies as ~ [R—R CTit y/ 4 and whose fastest growing 
wave number fc max ~ ^/const + (R — R cr it) near 
the instability threshold. Such an instability can 
also occur purely due to deposition dynamics ef- 
fects and is controlled by the atomic flux j at , 

• Type-II instability which is either driven by 
curvature-dependent surface-roughness processes 
or by synergistic interactions of preferential sput- 
tering and preferential curvature-dependent dis- 
placement/deposition. For the surface-roughness 
instability, the band of wavelengths and the fastest 
growing wavelength are characterized both by ~ 
\/R — i? C rit • Such type of an instability can also 
occur purely due to deposition dynamics effects. If 
both ion and deposition-driven effects act in induc- 
ing instability or one of them is negligible, both the 
unstable wavelength band and the fastest growing 
wavelength are given by ~ y/~R (or ~ \fja± ) . There- 
fore the wavelength diverges when approaching the 
instability threshold; 

• The oscillating or type-Ill instability occurs due 
to surface-roughness effects. It is characterized by 
a wavelength band width ~ (R — i? C rit) 1 ^ 4 and a 
finite fastest growing wave number ~ y/R— R cl n 
near the threshold. This instability can occur only 
in the presence of ion irradiation. 

Oppositely, ion irradiation and deposition can act in a 
mutually excluding manner resulting in the surface sta- 
bilization. 

In addition to causing surface instabilities, surface to- 
pography does not only influence the local redistribu- 
tion of the deposited material but also influences the net 
growth rate and the average steady state surface compo- 
sition. This implies that curvature effects must be con- 
sidered in the deposition balance. 

The reason for the fact that ion-induced instabili- 
ties can occur via either surface roughness processes or 
via composition-roughness feedback interactions is that 
we have considered only alloy systems. This translates 
into the term ~ d 2 <j)/dx 2 being multiplied by a posi- 
tive diffusion constant D* which acts in homogenizing 
the spatial alloy component distribution. In the case 
of spinodal decomposition, the coefficient D* becomes 
negative ; 10 ' 29 ' 77 and therefore the phase separation oc- 
curs spontaneously without any roughness-composition 
feedback interactions. 

In the case of nucleation and growth, one has to 
consider nucleation events as finite, not infinitesimal, 
fluctuations.— Therefore one has to deal with a non- 
linear theory from the beginning. This can be done by 
introducing the ion-induced and film growth terms of Eq. 
(fT8|) into the Cahn-Hilliard-Cook type equations.- 8 Both 
these cases (spinodal decomposition and nucleation and 
growth) are more complex and therefore requires more 
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sophisticated approaches. The present work provides a 
basic framework to address this problem by employing 
ion-to-atom arrival ratio R as a key control parameter to 
induce patterning during film growth. We believe that 
this parameter will play a key role for the phase separat- 
ing systems. 

There is an additional time dependence originating 
from the ion-induced surface drift which has not been 
considered in the analysis. This results in a pattern lat- 
eral shift as a function of film thickness. While the de- 
tailed analysis requires separate theoretical and experi- 
mental studies, we note that such drift effects combined 
with a possibility to rotate the sample during the film 
growth opens opportunities to sculpt complex 3D struc- 
tures such as chevrons or helices. Usually the growth 
of such structures is achieved by using glancing angle 
deposition . 79 i 80 Protection from the degradation of such 
structures requires a post-growth filling of the spaces be- 
tween the sculpted nano(micro) objects. The approach 
outlined in this work would inherently result in an en- 
capsulated nanostructure. 

In general, the terms describing other kinetic instabil- 
ities during film growth such as Ehrlich-Schwoebel insta- 
bility have similar mathematical form^ Thus they are 
also inherently described by the model presented in this 
manuscript. The physical interpretation is different and 
has been omitted in the discussion of the present pa- 
per. However, such instabilities are relevant in both the 
contexts - film growth 25 and ion-induced patterning 42 . 
Therefore, the reiteration of the general findings of this 
work in these contexts might also provide some new phys- 
ical insights on how deposition and ion-irradiation effects 
can synergystically affect the nanopattern formation. 

The results of this study present an alternative way to 
control surface patterns, especially if only a thin nanos- 



tructured layer of different material than that of a bulk 
substrate is needed. During the film growth the bulk usu- 
ally remains 'frozen'. Therefore, such a composition pat- 
tern formation on the surface during alloy growth trans- 
lates into an ordered 3D nanocomposite structure. De- 
pendent on the type of stationary or oscillating instabil- 
ity, this can result in laterally or even vertically ordered 
nanostructures. It should be noted that a self-organized 
multilayer structure is expected to occur neither due to 
sequential deposition nor due to sub-surface-driven phase 
separations^ but is purely an ion- induced surface effect. 

The effects analyzed here are of physical origin. 
The large energy delivered by ions provides efficient 
means to compete with other less energetic and material 
specific processes. We believe that if the ion irradiation 
is intensive enough (ion-to-atom arrival ratio R is 
above critical), such effects can be observed in many 
alloy systems. This enables a material-independent 
structural design approach to sculpt the morphology 
at the nanoscale. As different material systems exhibit 
different range of properties, this is of possible relevance 
for a wide range of applications. 
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FIG. 4: The phase diagram for preferential relocation to sputtering ratio of species A (as Fig. [3^,). The deposition effects are 
included considering the case of preferential curvature-dependent deposition to sputtering ratio of species A. This shifts all axes 
according to Eq. (|48|l by the corresponding Q values of Eq. (|49[1 . 



